Resonant activation driven by strongly non-Gaussian noises. 
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The constructive role of non-Gaussian random fluctuations is studied in the context of the passage 
over the dichotomously switching potential barrier. Our attention focuses on the interplay of the 
effects of independent sources of fluctuations: an additive stable noise representing non-equilibrium 
external random force acting on the system and a fluctuating barrier. In particular, the influence 
of the structure of stable noises on the mean escape time and on the phenomenon of resonant 
activation (RA) is investigated. By use of the numerical Monte Carlo method it is documented that 
the suitable choice of the barrier switching rate and random external fields may produce resonant 
phenomenon leading to the enhancement of the kinetics and the shortest, most efficient reaction 
time. 

PACS numbers: 02.50.-r, 05.10.-a, 05.90.+m, 82.20.-w 



X 



INTRODUCTION 



The fluctuating barrier problem seems to be of rele- 
vance especially for the kinetics of biophysical systems 
[E S S H whose time-evolution is governed by coupling 
to protein molecules undergoing constant conformational 
transitions. In particular, it is by now well established 
that membrane proteins that conduct many chemical and 
physical processes can exist in several conformational 
states whose distinct electrical properties may be respon- 
sible for the efficiency of ionic transport through biolog- 
ical channels. In such cases, the assumption of random 
fluctuations in the membrane potential has been shown 
to be sufficient to explain experimental data on 
active ionic flow and resonance between internal fluctu- 
ations and external driving AC field. The enhancement 
or optimisation of the action of the external regular field 
by stochastic fluctuations gained the term stochastic res- 
onance (SR) and became an issue of vivid experimental 
and theoretical studies [|| over the past 20 years. Al- 
though the SR phenomenon is commonly described as a 
cooperative effect in which small periodic influence en- 
trains external random noise, similar constructive effects 
of noises under nonequilibrium constraints can be ex- 
pected in "ratchet" systems JtJ where spatially uniform 
mean-zero symmetrical (random or periodic) time depen- 
dent forces interact with an underlying spatial anisotropy 
to induce motion. Another resonant- type, noise-induced 
behaviour [H IS 01 El 1S resonant activation (RA), the 
phenomenon in which an optimal fluctuation rate exists 
which minimizes the mean first passage time for the es- 
cape over the fluctuating barrier. The generic theoret- 
ical models analysing the aforementioned resonant phe- 
nomena are usually based on a Langevin equation ap- 
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proach assuming the overdamped limit 0, 0, El fl2| . 
Accordingly, the influence of the external thermal bath 
of the surroundings on a Brownian particle is described 
in such an equation by time-dependent random force 
which is commonly assumed to be represented by a white 
Gaussian noise. That postulate is compatible with the 
assumption of a short correlation time of fluctuations, 
much shorter than the time-scale of the macroscopic mo- 
tion and assumes that weak interactions with the bath 
lead to independent random variations of the parame- 
ter describing the motion. In more formal, mathematical 
terms Gaussianity of the state-variable fluctuations is a 
consequence of the Central Limit Theorem which states 
that normalized sum of independent and identically dis- 
tributed (i.i.d) random variables with finite variance con- 
verges to the Gaussian probability distribution. If, how- 
ever, after random collisions jump lengths are ruled by 
broad distributions leading to the divergence of the sec- 
ond moment, the statistics of the process changes sig- 
nificantly. The existence of the limiting distribution is 
then guaranteed by the generalized Levy-Gnedenko fl3| 
limit theorem. According to the latter, normalized sums 
of independent, identically distributed random variables 
with infinite variance converge in distribution to the Levy 
statistics. At the level of the Langevin equation, Levy 
noises are generalization of the Brownian motion and de- 
scribe results of strong collisions between the test par- 
ticle and the surrounding environment. In this sense, 
they lead to different models of the bath that go be- 
yond a standard "close-to-equilibrium" Gaussian descrip- 
tion 0HHG3. 

In this communication we present numerical results 
for the mean first passage time over a fluctuating bar- 
rier for the model system with a linear potential sub- 
ject to Markovian dichotomous fluctuations and additive 
Levy noises. We report here on a new phenomenon of 
reappearance of RA in a model system whose pattern 
of the most efficient kinetics has been smeared out and 
destroyed by the presence of driving a-stable noise. In 
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particular, we demonstrate that a congruent tuning of 
noise asymmetry with its stability index responsible for 
scaling properties of the noise increments can optimize 
the passage over the conformationally changing barrier. 



II. THE MODEL 

An overdamped Brownian particle is moving in a po- 
tential field between absorbing (x = 1) and reflecting 
(x = 0) boundaries, in the presence of noise that modu- 
lates the barrier height. Time evolution of a state vari- 
able x(t) is described in terms of the generalized Langevin 
equation 

^ = —V'(x) + grj(t) + V2((t) = -V' ± {x) + V2((t), (1) 

where prime means differentiation over x, is a white 
Levy process originating from the contact with non- 
equilibrated bath and ij(t) stands for a Markovian di- 
chotomous noise of intensity g taking one of two possible 
values ±1. Autocorrelation of the dichotomous noise is 
set to (fa(t) - (v))(v(t') - (r?))) - exp(-2 7 |i - t'\). For 
simplicity, throughout the paper a particle mass, a fric- 
tion coefficient and the Boltzman constant are all set to 
1. The time-dependent potential V±(x) is assumed to be 
linear with the barrier switching between two configura- 
tions with an average rate 7 

V ± (x)=H ± x, g= H -~ H+ . (2) 

Both £ and 77 noises are assumed to be statistically inde- 
pendent. 

The initial condition for Eq. JQ) is x(0) = 0, so that ini- 
tially particle is located at the reflecting boundary with 
equal choices of finding a potential barrier in any of two 
possible configurations (±). The quantity of interest is 
the mean first passage time (MFPT), the average time 
which particle spends in the system before it becomes 
absorbed. In the approach applied herein information 
on the MFPT is drawn from the statistics of numerically 
generated trajectories satisfying the generalized Langevin 
equation JTJ and the examined MFPT is estimated as a 
first moment of the distribution of first passage times 
(FPT) obtained from the ensemble of simulated realiza- 
tions of the stochastic process in question. More pre- 
cisely, frequency of FPT is calculated as frequency of the 
time-duration r of the events (trajectories) that have ini- 
tially started from x — and reached x — 1 for the first 
time. The MFPT is estimated as (r) over the frequency 
distribution. Despite the fact that in Eq. x itself is a 
random variable distributed according to some unknown 
stable distribution, the FPTs are distributed according 
to some unknown probability distribution possessing all 
moments. 



III. RESULTS OF SIMULATIONS 

Generally, for systems driven by Gaussian noises, 
noise-induced phenomena are studied solely in terms of 
tr 2 , the intensity of the white noise 0, Ug. Thermo- 
dynamic definition of this noise-source relates a 2 to the 
system temperature. It is not so, however, for dynamic 
systems perturbed by Levy stable noises, for which an 
appropriate choice of a control parameter is less ob- 
vious. It is caused by the fact that stable distribu- 
tions are characterized by a four-parameter family: the 
stability index a (a € (0,2]), skewness parameter (3 
(P G [-1,1]), the shift /J, {n e R) and scaling a (a e K + ) 
parameters, respectively (for the definition of Levy mea- 
sure, see Appendix A). The PDFs of stable distributions 
demonstrate asym ptotic p ower-law behaviour according 
to L(0 w Id" 1 "" USUI]. The parameter /it controls the 
position of the PDF modal value, whereas a, similarly 
to the Gaussian case, scales the width of the distribu- 
tion. At the first glance, the parameter a is somehow 
connected to the system temperature, although the re- 
lation is more subtle since stable additive noises can be 
expected in nonequilibrium situations where the defini- 
tion of temperature looses its common sense |0, 0] . 

For the purpose of simulations the value of fj, has been 
arbitrary set to and a = 1/ y/2. Such a choice of /1 and 
a reconstructs standard normal distribution iV(0, 1) for 
a = 2, therefore allowing for comparison of the results 
obtained in this study with the former ones Q, 0, 0] , 
where the Gaussian measures of the £ noise have been 
used. Remaining parameters (a and f3) have taken val- 
ues from the allowed range with a lower limit set up for 
the a value (a < .2 have not been investigated). The RA 
phenomenon has been examined for various barrier se- 
tups: The potential barrier has been switching between 
different heights H + and H- and as values of H± the 
sets ±8; (0, 8); and (4, 8) have been chosen. They cor- 
respond to the changing of the potential barrier between 
the barrier and the well (H± — ±8), the barrier and "no 
barrier" (H + = 8, H- = 0) situations and to two barriers 
with different heights (H+ = 8, H- = 4). 

In order to test the implemented numerical procedures 
(cf. Appendix A), MC solutions of Eq. for a = 2 have 
been compared with the exact solutions to the backward 
Fokker Planck Equation [H ITU]. The numerical results 
(see Fig. have shown a perfect agreement with the 
exact solutions. 

Subsequent figures demonstrate results of the analysis 
for driving stable noises with various values of the (a, (3) 
parameters. For the sake of clarity of the comparison, 
they have been grouped in three classes corresponding to 
different scenarios of barrier-alteration. The first group 
(Figs. [21 and |3) relates to the potential barrier switching 
between H± — ±8, Figs. 0] and describe numerical re- 
sults obtained for the barrier switching between H + = 8 
and H_ = and, eventually, Figs. HO and [7] refer to simu- 
lations performed for the potential barrier switching be- 
tween two different barriers, i.e. H+ = 8, H- = 4. In 
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FIG. 1: MFPT(7) for linear potential barriers switching 
between different heights H±: (+) H+ = 8, H- = 4; 
(x) -H+ = 8, H- = 0; (*) H+ = 8, H- = -8. The driv- 
ing noise is Levy stable noise with a = 2, i.e. Gaussian noise. 
Solid lines represent exact result constructed by direct inte- 
gration of the backward Fokker-Planck equation. Numerical 
results were obtained by use of Monte Carlo simulation of 
Eq. Q with time step dt = 10~ 4 and averaged over N = 10 3 
realizations. Error bars represent deviation from the mean 
and usually remain within the symbol size. 



each class, in the first (left) plot (cf. Figs. [21 EJandEJ) 

the results for fixed f3 {j3 = 0), i.e. for the symmetric 
stable noises are displayed while the right panel presents 
results for a fixed value of a (a = 0.5), thus demonstrat- 
ing the effect of the asymmetry of the driving noise. The 
second set of plots in each class (cf. Figs. El El and 0) 
presents results for the MFPT(a, 7) for a fixed f3 param- 
eter (/3 = 1) while a corresponding right panel displays 
sample cross-sections of the MFPT(a,7) surface. 

Closer examination of Figs. and fallows to conclude 
that a typical, non-monotonic character of the MFPT 
curve is preserved for all values of the a parameter. The 
only visible changes are registered in the asymptotic be- 
havior and in the depth of the MFPTs curves. In partic- 
ular, the values of MFPTs for the Gaussian case (a = 2) 
are higher than for other values of a thus demonstrat- 
ing that the heavier tails in PDFs (and therefore higher 
probability of experiencing larger fluctuations) of driving 
random forces facilitate the transport over the fluctuat- 
ing barrier. Furthermore, as it can be seen in the right 
panel of Fig. |5) for o = 0.5 the RA phenomenon is ob- 
served for all values /3. 

In Fig. |2| MFPT(a, 7) surfaces for asymmetric sta- 



ble noises (f3 = 1) are presented. Interestingly, the 
MFPT(a, 7) surface for asymmetric {[3 = 1) driving noise 
displays an unexpected behavior. Namely, a gradual 
change of the a parameter (from a = 0.2 to a = 1.1) 
results in suppression of the RA phenomenon that disap- 
pears for values of a ~ 1 but becomes further recovered 
for a > 1 approaching a representative for the Gaus- 
sian case. By comparison (cf. Fig. El left panel), in 
the cases when the dynamics flickers between the motion 
over the erected potential barrier and a free (symmet- 
ric) Levy flight, the characteristic resonant shape of the 
MFPT curve disappears with a decrease of the stability 
index a, i.e. noises with heavier asymptotic tails destroy 
the RA. Contrary to the H± = ±8 case, for a — 0.5, the 
RA phenomenon is not observed, independently of the 
value of P (cf. Fig. El right panel). However, for asym- 
metric stable noises (Fig. EJ the same kind of behavior 
which has been observed for H± = ±8 can be noticed. 
Here again, it is possible to detect an optimal a-value at 
which activation processes are more efficient (they lead 
to lower values of the MFPT than for all other choices 
of the a parameter). Unlike in the H± = ±8 case, the 
RA does not take now place for very small values of a 
but it reappears for a ~ 1 when the resonant shape of 
the MFPT curve becomes clearly visible. Afterwards, 
for moderate a the RA phenomenon disappears and, as 
expected, becomes again induced for a = 2. 

Figs. El and □ display results for H + = 8, H- = 4. As 
can be observed in the left panel of Fig. El the Gaussian 
behavior of the MFPT curve disappears with decrease of 
the stability index a for symmetric stable noises. The 
shape of the MFPT curve typical for the Gaussian case 
vanishes much faster than for the H+ = 8, H- = sce- 
nario. For a heavy tailed (a = 0.5) external noise driving 
(right panel of Fig. EJ) the RA is not visible. 

Analogously to the both previously studied barrier se- 
tups with asymmetric stable driving noises (Fig. 0) there 
exists an optimal a for which activation processes are 
most efficient. Also, the RA phenomenon is inhibited 
at low values of a, reappears for a ~ 1 and smears out 
for moderate a-s before its recovery becomes effective at 
a = 2. 

By inspection of Figs. El El and it can be concluded 
that the same kind of resonant behavior (described by 
a nonmonotonic shape of the MFPT curves) can be ob- 
served not only for equilibrium fluctuations but also for 
heavy-tailed noncquilibrium fluctuations modelled by an 
asymmetric noise with the stability index a ~ 1. Closer 
analysis of the MFPT (a, 7)-surface cross-sections shows 
that the phenomenon disappears for very small a (except 
of the H± = ±8 case), becomes induced for a w 1, disap- 
pears for a > 1 and renters again in the Gaussian limit. 
Note that the Gaussian cases (a = 2) have not been plot- 
ted in the figures with cross sections due to the fact that 
values of the MFPT for the Gaussian driving noise are 
significantly higher than for all other stable noises with 
smaller values of a. Therefore Fig. presents separately 
results of numerical analysis for a = 2. 
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FIG. 2: MFPT(a,7) for f3 = (left panel) and MFPT(/3, 7) for a = 0.5 (right panel) for the linear potential barrier switching 
between H± = ±8. The results were calculated by direct integration of Eq. £Q| with the time step dt = 10~ 4 and averaged over 
TV = 10 3 realizations. jf 
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FIG. 3: MFPT(a, 7) for H± = ±8 and (3 = 1 (left panel) and sample cross-sections MFPT(7) for various a: (+) a = 0.2; 
(x) a — 0.8; (*) a = 0.9; (□) a — 1.1 (right panel). The results were calculated by direct integration of Eq. 10 with the time 
step dt — 10 -4 and averaged over TV = 10 3 realizations. Error bars represent deviation from the mean. 



Examination of left panels in Figs. [21 El and allows to see how the typical Gaussian behavior of the MFPT 
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FIG. 4: MFPT(a,7) for /3 — (left panel) and MFPT(/3, 7) for a = 0.5 (right panel) for the linear potential barrier switching 
between H+ — 8, H- — 0. The results were calculated by direct integration of Eq. Q with the time step dt — 10~ 4 and 
averaged over N = 10 3 realizations. 

J* 
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FIG. 5: MFPT(a, 7) for H+ = 8, H- = and = 1 (left panel) and sample cross-sections MFPT(7) for various a: (+) a = 0.2; 
(x) a = 0.8; (*) a = 0.9; (□) a — 1.1 (right panel). The results were calculated by direct integration of Eq. Q with the time 
step dt — 10 -4 and averaged over TV = 10 3 realizations. Error bars represent deviation from the mean. 



curves changes with the change of the stability index a. For H± = ±8 only changes observed are visible in the 
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FIG. 6: MFPT(a,7) for /3 — (left panel) and MFPT(/3, 7) for a = 0.5 (right panel) for the linear potential barrier switching 
between H+ — 8, H- = 4. The results were calculated by direct integration of Eq. Q with the time step dt — 10~ 4 and 
averaged over N = 10 3 realizations. 
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FIG. 7: MFPT(a, 7) for H+ = 8, H- = 4 and /3 = 1 (left panel) and sample cross-sections MFPT(7) for various a: (+) a = 0.2; 
(x) a = 0.8; (*) a = 0.9; (□) a = 1.1 (right panel). The results were calculated by direct integration of Eq. Q with the time 
step dt — 10 -4 and averaged over TV = 10 3 realizations. Error bars represent deviation from the mean. 



asymptotic behavior of the MFPTs. For other cases un- der the study (H + = 8, H- = 0; H+ = 8, H- = 4), 
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the typical shape of the MFPT curve disappears with tial under the influence of periodic signal and a colored, 
decreasing a and the RA phenomenon vanishes. The non-Gaussian noise, 
process is the most rapid for H + = 8, if_ = 4 (cf. left 
panel of the Fig. EJ) while for H + — 8, H- = it slows 

down and becomes smoother (cf. left panel of the Fig.©. V. q-STABLE RANDOM VARIABLES AND 

a-STABLE LEVY MOTION PROCESSES 



IV. SUMMARY 

We have considered a thermally activated process that 
occurs in a system coupled to a non-Gaussian noise 
source introduced by a non-equilibrated thermal bath. 
Another external stochastic process is assumed to be re- 
sponsible for dichotomous fluctuations of the potential 
barrier which has been modeled by the linear function 
with a varying slope. 

The applied procedures enable to study the RA phe- 
nomena for all possible values of a, (3, a, /i mimicking sta- 
tistical properties of the driving stable noises. In partic- 
ular, for a = 2 (and any value of /?), the a-stable noise 
generated by the recipe given by Eq. is a Gaussian 
noise and numerically constructed results are in a perfect 
agreement with the exact results obtained by numerical 
integration of an appropriate backward Fokker-Planck 
equation [llj. In this case the RA takes place for all 
barrier setups under consideration. For a = 0.5, (3 — 1 
and a = l,/3 = formerly studied [l7| Levy-Smirnoff 
and Cauchy cases were recovered. Results of numerical 
analysis presented in this paper agree with previously 
constructed solutions |l7j . 

In the case of driving stable noises, the RA phe- 
nomenon is most resistant for the H± = ±8 scenario 
where it is visible for all noises under the study. For 
symmetric stable noise it is possible to observe how the 
Gaussian behavior of the MFPT curves changes and the 
RA phenomenon in the system switching between two 
barriers {H+ — 8, H— = 4) or barrier and no barrier 
= 8, H- = 0) vanishes. A novel behavior has been 
registered for nonsymmctric stable noises: For very small 
values of the stability index a, the RA is not effective 
(except of the H± = ±8 case). The phenomenon reap- 
pears by increasing a to « 1, vanishes for a > 1 and 
becomes again observable for a = 2 (the Gaussian case). 
Moreover, in the a space another kind of phenomenon, 
allowing to choose an optimal value of a for which the 
MFPTs are smaller than for other values of stability in- 
dex a, occurs. The fact that the shape of MFPT curve 
for very heavy-tailed, skewed distributions representing 
non-Gaussian noises arising from the contact with not 
equilibrated bath is similar to the MFPT curve in the 
presence of the Gaussian fluctuations arising from the 
contact with equilibrated thermal bath is very interest- 
ing and makes the problem of recognition of the underly- 
ing noise more complicated. This indicates that more de- 
tailed study of "resonant phenomena" in non-equilibrium 
situations might be required, as has been recently also 
argued in other studies [2(j relating to the detection of 
the stochastic resonance in a generic double well poten- 



The a-stable variables are random variables for which 
the sum of random variables is distributed according to 
the same distribution as each variable, i.e. 



all + bX 2 = cX + d, 



(3) 



where = denotes equality in a distribution sense. Real 
constants c, d in Eq. © allow for rescaling and shifting 
of the initial probability distribution. The characteris- 
tic function of the stable distribution can be parameter- 
ized in various ways. In the usually chosen L a ^{Q\<j, fj,) 
|!9ll2l| parameterization, a characteristic function of the 
Levy type variables is given by 



exp 



1 - z/3sign(/c) tan —J + i/zjfe 



i(3ka a tan , for a ^ 1, 



exp 



-a\k\ 1 + sign(fc) In Ifcl 

7T 



for a = 1, 



(4) 



with a € (0,2], f3 e [—1,1], c G (0,oo), /i G (-00,00) 
and 4>{k) defined in the Fourier space 



(5) 



The above parameterization Q is continuous, in the 
sense that 

7rrv 2 

lim a a (|fc| Q sign(fc) - fc) tan — = cr|/c|sign(fc) In |fc| 

a— >-l 2 7T 

(6) 

for every a and k. Analytical expres- 

sions for stable probability distributions 
L a ,p{Ct a i A*) are known only in few cases: for 
a = 0.5, (3—1 resulting distribution is Levy-Smirnoff 



x exp 



2(x - n) 

for a — 1, /3 = it results in the Cauchy distribution 



1 



(8) 



whereas for a = 2 with arbitrary (3 the PDF is Gaus- 
sian. Characteristic feature of distributions L Qj ^(^;cr, /i) 
is existence of moments of order a, i.e. the integral 
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J n ^ o i Q .^(C; cr, n)( a d( is finite. This statement results in 
the conclusion that the only stable distribution possess- 
ing the second moment is a Gaussian one; for all other 
values of a the variance of a stable distribution diverges, 
and for a < 1 also the average value does not exist. 

For the purpose of analysis, the corresponding 
Langevin equation Q has been simulated by use of the 
appropriate numerical methods. Position of the particle 
has been obtained by direct integration of Eq. (JTJ 



x(t) 



t p-t 
[V'(x(s))-gr)(s)]ds+ / dL a , fi (s) 

to J to 

t pt 
V±(x(s))ds + / dL a> p(s). 

to J to 



In general 0,0, El, the L a ,p measure in Eq. @ can 
be approximated by 



N-l 



to 



f(s)dL a , (s) « ^2f(iAs)M a , l3 {{iAs,(i + l)As)) 

N-l 

4 ^ /(iA s )A s V« ft , (io) 



i=0 

distributed 



with 
NAs = 



the PDF 

t — to and 



where q is 
La.pfcv = 7f' A 1 = 0), 
M a ^([iAs, (i 4- l)As)) is the measure of the inter- 
val \iAs, (i + I) As). 

Random variables ? corresponding to the characteristic 
function (lean be generated using the Janicki-Weron 
algorithm |l9| . For a ^ 1 their representation is 



D 



a>,{3,(T ~ 



sin (a(F + C a ,p)) 

(cOs(V))* 

cos(V" - a(V + C a ,f,)) 



w 

with constants B,C,D given by 
B a ,/3,cT,n = A 4 



-Ba,/3,CT,/i! (11) 



/3anan(^), 



(12) 



arctan (/3tan(^)) 



(13) 



D a ,0,a — °~ 



( ( n ,™. 

cos ^arctan ^ptan(-^-) 



(14) 



(9) For a = 1, ? can be obtained from the formula 



2(T 



<T = 



a, ft > 



(15) 



with 



H - 



-pVln(o-) 



(16) 



In the above equations V and W are independent ran- 
dom variables, such that V is uniformly distributed in 
the interval (— -| , ^) and W is exponentially distributed 
with a unit mean |l9j . Numerical integration has been 
performed for fi = with increments of AL a p sampled 
from the strictly stable distributions |2l|, i.e. in gen- 
eral stable distributions generated via Eqs. (|ll|) - (|16[l the 
value of fi is shifted such that effectively B a ^^^ is equal 
to 0. Due to documented |21j numerical instability of all 
calculations with a-stable random variables for the set 
(a = 1,|3 ^ 0), this case has been excluded from the 
analysis. Algorithms as described above have been posi- 
tively tested for independence on the choice of the inte- 
gration step (the evaluation of stochastic integrals with 
time steps dt — 10 -4 and dt = 10~ 5 yielded the same 
results) . 
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